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We introduce a two-flavor quark-meson-diquark model for two-color QCD and its extensions to 
include gauge-field dynamics as described by the Polyakov loop. Grand potential and phase structure 
are being studied both in mean-field approximation and with the functional renormalization group. 
The model provides an explicit example for the importance of baryonic degrees of freedom: When 
they are omitted, the phase diagram closely resembles that of the corresponding (Polyakov)-quark- 
meson models for QCD, in particular including their critical endpoint. In order to reproduce the 
well established main features based on the symmetries and breaking patterns of two-color QCD, 
however, they must be included and there is no critical endpoint. The competing dynamics of 
collective mesonic and baryonic fluctuations is well described by the functional renormalization group 
equation in lowest order derivative expansion for the effective potential which we solve numerically 
on a two-dimensional grid in field space. 
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I. INTRODUCTION 

The phase diagram of Quantum Chromodynamics 
(QCD) is subject to enormous international research 
campaigns [HE]. In order to understand its main char- 
acteristic features such as the different phases of strongly 
interacting matter, the nature of the transitions between 
them, the existence and locations of critical points or 
perhaps approximate triple points and even multicritical 
points, it has proven to be very useful to deform QCD by 
not only varying the individual quarks' masses but also 
the numbers of their different flavors and colors. An im- 
portant example is the limit of infinitely many colors N c 
which inspired many qualitative descriptions of the QCD 
phase diagram [3]-[5] . One interesting aspect of this limit 
is that the baryon density becomes an order parameter 
for N c — > oo, in particular, also when the number of fla- 
vors Nf grows along with N c , i.e. for Nf/N c held fixed. 
In this paper we study two-color QCD, which shares this 
aspect of the large- N limit, here with iV c = Nf = 2. 
If one accepts that the phases of many-color QCD with 
N c ~ Nf — > oo have a bearing on the real world, it might 
therefore not be absurd, with due appreciation of all dif- 
ferences, to consider N c = Nf = 2 either. 

Quantum Chromodynamics with two colors (QC 2 D) 
has been well studied for many years within chiral ef- 
fective field theory and random matrix theory [51414] . in 
lattice simulations T5H22] , and the Nambu-Jona-Lasinio 
model In this paper we formulate a Polyakov- 

quark-meson-diquark (PQMD) Model for studying the 
phase diagram of QC2D with the functional renormal- 
ization group, including fluctuations due to collective ex- 
citations. 
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The most important differences between two and three 
colors follow from the special property of the SU(2) gauge 
group of QC2D: Its representations are either pseudo-real 
or real which leads to an antiunitary symmetry in the 
Dirac operator [7]. As a result, the fermion determinant 
remains real for non- vanishing baryon or quark chemical 
potential, \i ^ 0, as it does for adjoint quarks in any- 
color QCD, or in the G2 gauge theory with fundamen- 
tal fermions also, for example. Thus, for an even num- 
ber of degenerate fundamental quark flavors in QC2D 
there is no fermion- sign 'problem and the phase diagram is 
amenable to lattice Monte-Carlo simulations. Symmetry 
considerations, lattice simulations and non-perturbative 
functional continuum methods together should therefore 
allow us to understand the phase diagram of this theory 
completely. A combined effort towards this goal will be 
very worthwhile in particular because it will help to bring 
the functional continuum methods to a level at which 
they can reliably be applied, with the necessary adjust- 
ments, also to real QCD where lattice simulations suffer 
from the infamous fermion-sign problem [33 . 

Another consequence of the pseudo-reality is the Pauli- 
Giirsey symmetry which allows to combine quarks and 
charge-conjugated antiquarks into enlarged flavor multi- 
plets. As a result, for vanishing chemical potential and 
quark mass, fi = m q = 0, the usual SU(Nf) x SU(Nf) x 
U(1)b chiral and baryon number symmetries are replaced 
by an extended SU(2Nf) flavor symmetry which is (spon- 
taneously) broken by a (dynamical) Dirac mass down to 
the (2Nf + l)Nf dimensional compact symplectic group 
Sp(Nf). For Nf — 2 the extended flavor symmetry 
group SU(4) and its Sp{2) subgroup are locally isomor- 
phic to the rotation groups SO(6) and SO(5), respec- 
tively. The coset is given by S 5 , the unit sphere in six 
dimensions, and a spontaneously generated Dirac mass 
will lead to five Goldstone bosons, the three pions plus a 
scalar diquark-antidiquark pair. 

Moreover, for N c — 2 these color-singlet scalar di- 
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quarks play a dual role as bosonic baryons at the same 
time. While this thus represents the perhaps most im- 
portant difference as compared to the real world, it also 
makes it much easier to investigate the effects of baryonic 
degrees of freedom on the phase diagram in functional 
approaches. In that sense our model can be considered 
as a first step towards their inclusion in a 'quark-meson- 
baryon' model for real QCD with three colors. 

For the same reason our model of QC2D provides a rel- 
ativistic analogue of the BEC-BCS crossover in ultracold 
fermionic quantum-gases, which has also been described 
successfully with functional renormalization group meth- 
ods 34. 35j. In contrast to non-relativistic models of the 
BEC-BCS crossover, an interesting additional constraint 
thereby arises from the Silver Blaze property |36j : When 
a relativistic chemical potential /1 is coupled to degrees 
of freedom with a mass gap A, at zero temperature, the 
partition function and hence thermodynamic observables 
must actually remain independent of the chemical poten- 
tial as long as fi < A. We will see that it is not trivial, in 
general, to implement this constraint in non-perturbative 
functional renormalization group studies, and that it can 
provide valuable extra information to devise intelligent 
truncations. 

Our main interest here, however, is to explicitly 
demonstrate the impact of baryonic degrees of freedom 
on the phase diagram by comparing the purely mesonic 
model, representative of typical three-color QCD model 
calculations, to the full quark-meson-diquark model. 

For this comparison we argue that it is more appropri- 
ate to think of the vacuum diquark mass as the baryon 
mass tub rather than the pion mass m,. In QC2D with 
its extended flavor symmetry they are the same, but the 
essential aspect of this assignment is that a continuous 
phase transition at zero temperature occurs at a critical 
quark chemical potential fi c — m,B/N c . Except for the 
scale separation between m n and in the real world, 
this transition can then be thought to correspond to the 
liquid-gas transition of nuclear matter in QCD with three 
colors which is of first order, involves the binding energy, 
and thus occurs somewhat below /1 = ms/iV c . 

As temperature increases the liquid gas transition 
ends, turning into a crossover with continuously vary- 
ing but nevertheless probably still relatively abruptly in- 
creasing baryon density along some narrow region. This 
rapid increase is generally expected to lead to the strong 
chemical-potential dependence of the chemical freeze-out 
line observed in heavy ion collisions at center-of-mass en- 
ergies below about 10 GeV per nucleon pair, the baryonic 
freeze-out [37l |38] . One might conclude that the phase 
transition line for diquark condensation, where a rapidly 
increasing baryon density spontaneously develops, would 
be the origin of a corresponding baryonic freeze-out line 
in two-color QCD, with N c = Nf = 2 arguably not nec- 
essarily further from reality than the large N c limits. As 
in the latter, one might then even identify a two-color 
version of quarkyonic matter [H [5TJ [35] . 

Finally, we would like to point out that our model, 



the functional renormalization group equations and the 
techniques to solve them have a broad scope of appli- 
cations beyond two-color QCD. One example is QCD 
with two light flavors at finite isospin chemical potential, 
which has been studied with the NJL model in mean- 
field plus random phase approximation (RPA) [39l 00] . 
There is a precise equivalence between the correspond- 
ing quark-meson model with isospin chemical potential 
and our quark-meson-diquark model of two-color QCD. 
Besides changing N c this merely involves reducing the 
number of would-be Goldstone bosons from five to three 
again, retaining only one of our degenerate pions, and 
reinterpreting the diquarks as the charged pions with 
isospin chemical potential [41 j . Similar models are also 
studied in the context of color superconductivity [12TFH] . 
The capacity to numerically solve functional renormaliza- 
tion group equations on higher dimensional grids in field 
space is generally useful for competing symmetries, as in 
a quark-meson model study of the axial anomaly with 
scale dependent 't Hooft couplings, for example |4"5] , 

The outline of this paper is as follows: In Sec. [E] we re- 
view the general features of QC2D such as its enlarged fla- 
vor symmetry and the possible symmetry breaking pat- 
terns in some more detail. Based on these symmetry 
considerations we then construct the Polyakov-loop ex- 
tended quark-meson-diquark Lagrangian for QC2D. In 
the next section, Sec. |III| we derive the thermodynamic 
potential of the (P)QMD model in mean-field approxima- 
tion, discuss so-called vacuum contributions, the Silver 
Blaze property and the relevance of pole versus screening 
masses for mesons and diquarks. The functional renor- 
malization group flow equations for the effective poten- 
tial in leading-order derivative expansion are derived in 
Sec. |IV| In this section we also calculate critical expo- 
nents which are consistent with the expected symmetry 
breaking pattern, investigate in how far mean-field re- 
sults are modified by fluctuations, and give a transpar- 
ent illustration of the importance of baryonic degrees of 
freedom for the phase diagram. As a byproduct we note 
that starting from a tricritical point, a region of first- 
order transition limiting the diquark condensation phase 
at larger chemical potentials as predicted from chiral per- 
turbation theory at next-to- leading order jjj], is also ob- 
served in the QMD model at the mean-field level. This 
first-order transition turns out to be a mean-field artifact, 
however. It is washed out by the fluctuations, and there is 
no sign of a tricritical point left, once the thermodynamic 
potential is obtained from its functional renormalization 
group flow. We draw our conclusions and present an out- 
look in Sec. [Vj Technical details can be found in several 
appendices. 



II. FLAVOR SYMMETRIES IN QC 2 D 

We begin this section with a short review of the ex- 
tended flavor symmetries of QC2D due to its Pauli- 
Giirsey symmetry, and the associated symmetry breaking 



3 



patterns. We then discuss a qualitative phase diagram for 
two-flavor QC2D and construct the quark-meson-diquark 
(QMD) model by a suitable vector coupling of quark bi- 
linear s to meson and diquark fields. 

A. Extended flavor symmetries and symmetry 
breaking patterns 

As all half-odd integer representations of SU(2), its 
fundamental representation is pseudoreal, which means 
that it is isomorphic to its complex conjugate repre- 
sentation with the isometry given by S — i(72, S 2 — 
— 10 Therefore, charge conjugation of the gauge fields in 
QC2D can be undone by the constant SU (2) gauge trans- 
formation S — i<T2- From now on we will use T a = a a /2 
for the color generators, with 

T aT = T a* = _ ST ag-l ^ 

and reserve cr, (t<) for the Pauli matrices in spinor (fla- 
vor) space. Together with the charge conjugation matrix 
C in spinor space, likewise with C 2 = — 1, and complex 
conjugation denoted by K one then defines an antiuni- 
tary symmetry T = SCK with T 2 = +1 (in a real color 
representation with S 2 = 4-1, one has T 2 = — 1, corre- 
spondingly) . This leads to the classification of the Dirac 
operator by the Dyson index j3 of random matrix theory 
O [7] , with P = 1 for fermions in the pseudoreal funda- 
mental color representation of QC2D (or f3 = 4 in the 
real color representations of SU(N)/Zn or G^)- 

Following [7J, we start from the standard kinetic part 
of the Euclidean QC2D Lagrangian, in the chiral basis, 

£kin = ■tplpip = ^lia^D^L - ip^ia^D^ipH , (2) 

with hermitian 7-matrices, = (— i, a), and ipR/L> ^ril 
as independent Grassmann variables with ^r/l = ^rJl- 
The covariant derivative is — + iAn, and the cou- 
pling is absorbed in the gauge fields = A®T a . 

The two terms in ^ get interchanged under the anti- 
unitary symmetry T. If we apply it only to the second 
term, say, by using (— i<r 2 ) for the chiral i?-component of 
the charge conjugation matrix C, i.e., changing variables 
to tpR = —i<J2Sip R and ip* R = —\a2S1pR, we can therefore 
reexpress 

Adn - Vha^D^ (3) 

in terms of the 2Nf 4-dimensional spinors "J" = (ipzi "4>r) T 
and = (ip^ L , V'j?)- Because it is now block diagonal, the 
SU(2Nf) symmetry in the space combining flavor and 



1 The irreducible representations of the proper rotations are real 
which means their complex conjugates are obtained from isome- 
tries S with S 2 = +1, just as those of the adjoint groups 
SU(N)/Zn or most of the exceptional Lie groups such as Gi- 



transformed chiral components is manifest in this form. 
With the same transformation of variables the quarks' 
Dirac-mass term becomes 

fni)i> = y(* T icr 2 5S * - y* T ia 2 SZ y*) , (4) 
where the symplectic matrix 

acts in the 2iV/-dimensional extended flavor space. 
An explicit (dynamical) Dirac mass therefore explic- 
itly(spontaneously) breaks the original SU(2Nf) down 
to the compact symplectic group Sp(Nf), sometimes also 
referred to as USp(2Nf) reflecting the fact that it is the 
intersection of the unitary U(2Nf) and the symplectic 
Sp(2Nf,C), the invariance group of So as bilinear form 
on complex 27V/-vectors. 

For Nf — 2 flavors the enlarged flavor symmetry group 
of QC2D is SU(4), not U(4) because of the axial anomaly, 
it replaces the usual chiral and baryon number symme- 
tries SU(2)l x SU(2)r x U(1)b- Just as this extended 
flavor SU(4) shares its 15 dimensional Lie algebra with 
the group of rotations in 6 dimensions, SO(6), its Sp{2) 
subgroup leaving the Dirac-mass term invariant has the 
10 dimensional Lie algebra of SO(5) (in fact they are both 
the universal covers of the respective rotation groups). 

Our brief review of the QC2D symmetries so far holds 
for vanishing chemical potential. For fi 7^ but m = 0, 
the SU(2Nf) symmetry is broken explicitly by fii/jjotp to 
SU(N f ) L x SU(N f ) R x U{1). This is also easy to see, 
from Eqs. ([2]), ([3]), as it amounts to introducing the term 
/x^7o^ = n^Bo^ with [7] 

B a = -7oS„ = _°J . (6) 

For Nf — 2, in terms of the rotation groups, this sym- 
metry breaking pattern is locally the same as SO (6) — >• 
50(4) x SO(2). 

When both /1 and m are non-zero, the unbroken flavor 
symmetry is of course given by the common subgroup 
SU(2)v x f7 (1) of the two limiting cases p, — >• 0, m ^ 
Oorm— ^0,^7^0 discussed above. Whether, as an 
approximate symmetry, it is more like the Sp(2) ~ 5*0(5) 
or like the SU{2) L x SU{2) R x U(l) ~ SO (4) x SO(2), 
naturally depends on the relative sizes of Dirac mass m 
and quark chemical potential \i. 

More precisely, it is an exact result of chiral effective 
field theory [7J, that for baryon chemical potential 
fiB = 2fi < m n the approximate chiral symmetry break- 
ing pattern remains that of \i = and the vacuum align- 
ment is (qq)-like with an approximate Sp(2) ~ SO(5) 
if m is sufficiently small, while (at zero temperature) 
for /is = 2\i > a diquark condensate develops and 
the vacuum alignment starts rotating from being (qq)- 
like to becoming more and more (gg)-like as /i is fur- 
ther increased. The chiral condensate then rapidly de- 
creases, chiral symmetry appears to get restored but the 
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increasing chemical potential reduces it to the approxi- 
mate SU(2) L x SU(2) R ~ 50(4) again. It is not the 
full flavor symmetry of the /i ^ 0, m — case discussed 
above because we have entered the diquark-condensation 
phase with spontaneous baryon-number breaking, corre- 
sponding to superfluidity of the bosonic baryons. 

Another exact result is that, at zero temperature and 
for )1b < Wtt, the onset of baryon condensation, the 
baryon density remains zero and the thermodynamic ob- 
servables must be independent of /i. Because this is far 
from obvious to verify explicitly in actual calculations, it 
has been named the Silver Blaze Problem [36) . In order 
to be able to excite any states at zero temperature, and 
with a gap in the spectrum, the relativistic chemical po- 
tential needs to be increased beyond the mass gap in the 
correlations to which it couples. Here, with a continuous 
zero-temperature transition at \xb = ra-K this gap is sim- 
ply given by the baryon mass in vacuum which because 
of the extended flavor symmetry in QC2D coincides with 
the pion mass, ms = m^. This latter property is of 
course special to N c = 2. The Silver Blaze property will 
hold, as it does here, up to a quark chemical potential 
of the order of m B /N c (reduced by 1/N C of the binding 
energy per nucleon when the transition is of first order), 
in general, however. 
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FIG. 1. Schematic phase diagram for QC2D in the parameter 
space of temperature T, quark mass m, and baryon chemical 
potential fin- 

At finite temperature, a qualitative picture emerges for 
the phase diagram as sketched in Fig. [l] The solid line 
in the T = plane there represents the continuous zero- 
temperature transition with diquark condensation which 
is of mean-field type. Because the quark mass m q scales 
quadratically with the pion mass, it will occur along a 
parabola m q oc p? B . The thick dashed lines represent the 
corresponding second-order transitions at finite temper- 
ature in fixed m q planes of the 0(2) universality. The 
thick line along the temperature axis is the magnetic 
first-order transition in the /is = plane which prob- 
ably ends in a multi-critical point. When viewed in the 
= plane, this is the critical endpoint in the 0(6) 
universality class for the chiral phase transition in QC2D 
with its extended SU (4) flavor symmetry. In the m q = 
plane, the vacuum alignment will always be (g^-like, for 
no-matter-how-small /1 > 0. Therefore, in this plane one 



only has the second-order 0(2) line which might end in 
the same point making it multi-critical. 

B. Quark- meson-diquark model Lagrangian 

The starting point of our model construction is the fla- 
vor structure of the standard chiral condensate and the 
quark mass term which is of the form \E' T Eo^'- It there- 
fore transforms under the full flavor SU (4) according to 
the six-dimensional antisymmetric representation in the 
decomposition 4 ® 4 = 10 ffi 6. 

The other components belonging to the same multiplet 
are obtained from transformations 

-> U^ , U = exp(i(9 a A a ) G SU(4)/Sp(2) . (7) 

Then, \ff T £ \E r -> v^E*, where, from Cartan's immer- 
sion theorem, the whole coset SU(A)/Sp(2) = S 5 is 
obtained in this way via E = U T Y, a U. The coset el- 
ements E are in turn parametrized by six-dimensional 
unit vectors n as E = riE, with E]Ej + E^E^ = 2Sij 

and E = (£ ,i£ AT a ) such that X a , a = 1 ... 5, form a 
basis for the coset generators [TT]. Thus, one verifies ex- 
plicitly that the vector \E ,T £<J' transforms as a (complex) 
six-dimensional vector under 5*0(6). 

A locally SU(2) C invariant linear sigma model La- 
grangian can therefore be defined by coupling the real 
50(6) vector of quark bilinears (4 ,T E^I'+h.c.) to a vector 
of mesonic fields <f) = (er, 7?, Re A, Im A) T formed by the 
scalar a meson, the pseudoscalar pions 7? and the scalar 
diquark-antidiquark pair A. This yields the Lagrangian 
(now including color and spinor components again), 

C a =¥ia^D^ + ^(* T iCT 2 5E* - ^* T icr 2 5E**)^ 

1 . . (8) 

where V(<fi) is the meson and diquark potential whose 
precise form will be specified later. A non-vanishing 
chemical potential couples not only to the quarks but also 
to the bosonic diquarks. Rewriting Eq. (j8| in terms of 
the original variables we obtain the quark-meson-diquark 
(QMD) model Lagrangian 

£qmd =i> + g{a + i7 5 7?f) - ^7°) ijj 

+ I (A*0/> t C 7 5 t 2 5V>) + A(V> t 7 5 T 2 5^*)) 

+ i((cV-2^)A)(cV + 2 M <>°)A* , ( 9 ) 

with C = 7 2 7° and a flavor- and color-blind Yukawa 
coupling g. With 

V($) = ±($ 2 -vY-ca, (10) 
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one obtains the corresponding 0(6) linear sigma model; 
and in the limit A — > oo, the bosonic part of £qmd 
is equivalent to the leading-order %PT Lagrangian of 
Refs. [7] . To see this explicitly it is best to start from the 
latter, use the explicit coset parametrization of as 
given above, and make the identifications v = f w = 2F 
and c = /tttoJ = 2Fm 2 r . It maybe worthwhile men- 
tioning that the coefficient of the leading term in [i of 
the xPT Lagrangian, (i 2 tr(Y,B T ^ B) with B = UB U\ 
which was fixed from gauging the flavor SU(A) in [6], 
here simply follows from — 2/i 2 |A| 2 as part of the kinetic 
term of the complex scalar diquark field A with chemical 
potential \xb — 2/i. This implies in particular, that the 
meson/diquark potential V{4>) itself, up to the explicit 
breaking by — cer, which needs to be only £0(4) x SO{2) 
invariant in general at finite fj,, must remain SO(6) invari- 
ant, however, at this leading order, C(/i 2 ), and therefore 
at 0(4> 2 ) in the fields, likewise. We can thus only have 
an 5*0(6) invariant mass term in V(4>). 

In the following it will be more convenient to rewrite 
the Lagrangian in terms of Nambu-Gorkov-like spinors 

* = ( rti> c ) ' wnere Vv O^g) denote the red (green) color 
components of ip and ip c = Cip T as in [55] . This yields 



£ 



QMD 



f 



1 



{d^f + v(<p) 



-((3 M -2^°)A)(9 M + 2^)A*, 



where 



C-l _ ( 0+£K<T+i7 5 ff"r)-7V 
°° ~ \ "97 5 A* 



97 A 
I +g ( a — i7 5 if t) +7° 



(12) 



Gauge field dynamics and confinement effects can be 
modeled also in QC2D by including a constant Polyakov- 
loop variable as a background field as in the NJL model 
|28j . and analogous to what is commonly done in the 
so-called Polyakov-loop-extended quark-meson models of 
three-color QCD [46-48 . To this end one introduces a 
constant temporal background gauge field = AqS^o 
which is furthermore assumed to be in the Cartan sub- 
algebra as in the Polyakov gauge, i.e., for SU(2) C simply 
given by Aq = T 3 2ao- This leads to the Polyakov loop 
variable 



$ = ^Tr c e ifiA ° = cos(/3a ), 



(13) 



to model a thermal expectation value of the color-traced 
Polyakov loop at an inverse temperature ft — 1/T, as 
an order parameter for the deconfinement transition at 
vanishing chemical potential. The covariant derivative 
d u — ififioAo leads to an additional contribution 



of the form 



itpj°T 3 2aoip which can be rewritten as 
— i^7°ao v I / in terms of the spinor field *Sf. Finally, we 
then arrive at the Polyakov-loop-extended quark-meson- 
diquark model (PQMD) Lagrangian, 



C 



PQMD 



= £ 



QMD 



i* 



7 U ao 
7°ao 



*+W Pol ($), (14) 



with £qmd defined in Eq. (11) and Up i($) is the 
Polyakov-loop potential [28] which is commonly fitted 
to lattice results, but which can also be computed with 
functional methods [HU [50]. In contrast to the three- 
color case the Polyakov-loop potential is a function of 
one single real variable $ here, even in the presence of a 
diquark condensate. 



III. MEAN-FIELD THERMODYNAMICS 

The grand potential in mean-field approximation is ob- 
tained by integrating over the quark fields and neglecting 
bosonic fluctuations. This means that all mesonic and 
diquark fields are replaced by their constant expectation 
values a = (<r), A = (A), A* = (A*) and 7? = (7?) = 0. 
In momentum space we then obtain, 



/•MF 
^-PQMD 

where 



= * 



(s- 1 

^0,MF 



i7°a 



°0,MF 



— ijS— 7 ft+gcr 
-97 6 A* 



* + V MF (<J, d 2 ) + W Po l($), 
(15) 

(16) 



97 A 



and V MF (<T,d 2 ) = (A/4) (a 2 + d 2 - v 2 ) -ca-2fi 2 d 2 with 
d 2 = |A| 2 is the bosonic effective potential. The last 



(11) term comes from the kinetic diquark part of Eq. (fill) 



and is included in the effective potential here. This term 
and the explicit chiral symmetry breaking by —ca break 
the SU(4) symmetry of the effective potential Vmf- The 
details of the parameter fixing and the values used in the 
numerical calculations are given in Appendix |A") 

The fermion-loop integration then yields for the grand 
potential £1, 



fi(T,/x) 



T J2 /(2^ Trh ( 5 0,MF 

V M f{v, d 2 )+U Pol (<i>), 



(17) 



where the trace runs over internal indices (Dirac-, flavor- 
and Nambu-Gorkov space) and we sum over antiperi- 
odic Matsubara modes v n — (2n + l)nT. The four dis- 
tinct eigenvalues of 7°S'(7mf are given by ±i?+ — 'w n and 



i-Bp — \v n with 



E* = ^g 2 d 2 



and 



(18) 



g 2 a 2 . 



The Matsubara sum can be performed analytically with 
the result 



tt(M 2 ,$) 



d*pr. 
2tt) 3 V 

2<$>e- l3E t + 
V MF (a,d 2 )+U Pol ($). 



(2- 
Tln ( 1 



e; + e- 



(19) 



')} 
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When the bosonic potential Vmf is replaced by M 2 (a 2 + 
d 2 ) - cer, with M 2 = g 2 /(4:G) and c = 2gm /(4G), this 
coincides with the Hubbard-Stratonovich transformed 
PNJL model result [25] with four-quark coupling G and 
current-mass parameter mo. Note that the model inde- 
pendent ~2/jL 2 d 2 term from chiral effective field theory, 
which is included in the bosonic part of the (P)QMD 
model, does not explicitly show up in the grand potential 
of the (P)NJL model. Minimizing the thermodynamic 
potential with respect to the constant mean fields tr, <i, $ 
leads to the gap equations, 



dVt _ on 
da dd 



an 



, 



(20) 



whose simultaneous solution yields the temperature and 
chemical potential dependent condensates a, d and <F 



A. Vacuum contributions 

The fcrmion-loop contribution to the grand potential 
in mean-field approximation, Eq. (19 1, contains an ul- 
traviolet divergent vacuum part. 



In the standard no- 
sea mean-field approximation one usually dismisses this 
vacuum contribution to the bulk thermodynamics. For 
some phenomenological consequences of this additional 
approximation and its influence on mean-field results, see 
Ref. [5TJ [55] and the references therein. Here we add 
an observation concerning this mean-field ambiguity of 
the quark-meson model when viewed as the d — > limit 
of the quark-meson-diquark model grand potential. Be- 
cause the Polyakov-loop contributions are irrelevant here 
we set $ = 1. 

For d — Eq. ( 19 1 superficially appears to reduce to 
the conventional mean-field expression for the grand po- 
tential of the quark-meson model |53j up to an overall N c 
in front of the fcrmion-loop contribution O g , 



n™-4N c T 

q 



d J p 

(2tt) 3 



i 



-0(e p ±M) 



(21) 



We illustrate the effect of the vacuum contribution for 
N c = 2 colors but without diquark condensate (i.e. with 
d = 0) in Fig. [2| To regularize the vacuum term f^ ac we 
employ a simple three-momentum cutoff A and assess the 
dependence of the phase structure on A. The parameters 
are fixed to reproduce an i/Vc-scaled f n — 76 MeV and 
a pion screening mass of m n — 138 MeV. In each case the 
sigma meson mass is adjusted so as to yield a common 
value for a chiral transition at [i = of T c sa 183 MeV. 
A = corresponds to the no-sea approximation. The de- 
pendence of the position of the critical endpoint (CEP) 
at fi c on the cutoff A is clearly visible in Fig. [2j With 
increasing A its location gets shifted towards larger chem- 
ical potentials and approaches the dimensionally regular- 
ized result [H] when A//i c is sufficiently large. 

More carefully, however, one observes that the fermion- 
loop in the no-sea approximation (f2™ c = when A = 0) 



A = 0MeV 




A = 500 MeV 

"•^^ A = 650 MeV 

^ s \ i ^ dim. reg. " 











0.5 1 1.5 2 2.5 

FIG. 2. Standard N c — 2 QM phase diagram: dependence 
of the location of the CEP on the vacuum-term cutoff A in 
comparison to dimensional regularization. 



in Eq. ( 21 ) does not tend to zero for T — > when /i > ga, 
but still contains temperature independent contributions 
from momenta with p 2 < [i 2 — g 2 o 2 . 

On the other hand, the d — > limit of Eq. ( 19 ) with 

\e p - M | (22) 



e: 



yields a grand potential of the quark-meson model for 
two colors which depends only on the chiral condensate a 
but which differs from the conventional expression by the 
appearance of the modulus of the quasi-particle energies, 



O(a) =-4 J -^{ep + fi+lep-fil 
+ 2Tln (l+ e ~^ + ^ 

2Tln (l+ e^l^l)} + V MF (^0) 



(23) 



The last two terms herein deserve to be called thermal 
now, as they do vanish at zero temperature for all /i. 
When vacuum and the thermal contributions are regular- 
ized in the same way we can recover the usual expression 
by means of the identity 

|x| +21n(l + exp(-|a;|)) = 21ncosh(x/2) + 21n2 . (24) 

This is for example the case in the NJL model if one 
chooses to regulate both thermal and vacuum parts with 
a three-momentum cutoff, but this is not what is usually 
done in the quark-meson model where the ultraviolet fi- 
nite thermal contributions are meant to be fully retained. 
The cutoff in the phase diagrams of Fig. [2] was applied 
only to 17™ c in Eq. (21), likewise. Otherwise the picture 
would change yet again. This is a bit of a grain of salt 
for the no-sea mean-field approximation in quark-meson 
models which is best motivated phenomcnologically as 
modelling the restoration of chiral symmetry at T = 
for large chemical potentials. Luckily, the problem is 
irrelevant altogether, once fluctuations are included via 
the functional renormalization group for which the quark- 
meson model shows its true uses. 
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Meanwhile, for the mean-field analysis of our quark- 
meson-diquark model including the possibility of diquark 
condensation with d ^ 0, we really have no option other 
than the perhaps anyway more natural splitting of ther- 
mal and vacuum contributions based on the modulus of 
the quasi-particle energy as the d 
QMD model mean-field grand potential, Eq. (j 1 9 ) 



our quark-meson-diquark model mean-field result with 
vacuum contribution from Eqs. (19 1, (20 1. 



limit of E in the 



If one considers the difference between cutting off the 
vacuum term in Eq. (21) as compared to the one in 



Eq. ( 23 ) as measure for the reliability of the calculation, 



one is led to conclude that the cutoff A in fi™ c should 
always be larger than the chemical potential fi. 

In the following we will continue to regulate vacuum 
terms with a sharp momentum cutoff mainly because di- 
mensional regularization, as applied to the three-color 
PQM model in [ST], due to the structure of E^, gets too 
complicated for a semi-analytic treatment here, with full 
diquark mean fields from the grand potential in Eq. ( 19 1. 



B. Diquark condensation 

Independent of the discussion in the previous section 
and of the influence of fluctuations, we know for two-color 
QCD that the quark-meson-model-like phase diagrams of 
the form as those in Fig. [2] are wrong. The exact chiral 
effective field theory results [H [Jj from the symmetries 
and breaking patterns as reviewed in Sec. |II A| tell us 
that we must include the diquark condensate along with 
the chiral condensate and base our mean-field analysis on 
Eqs. ( 19 1 , ( 20 ) in order to describe the superfluid diquark 
phase starting at a critical line /i c (T) with (U c (0) = 
(or fi B = m B )- 




FIG. 3. Condensates at T — (NJL parameter values from 
|26| ; lattice data from |T7j ; linear sigma model with = 138 
MeV and m a = 680 MeV). 

The resulting chemical potential dependence of the chi- 
ral and diquark condensates at zero temperature is shown 
in Fig. [3] where we compare the prediction from leading 
order chiral perturbation theory [7J , the NJL j^Hj and the 
linear sigma model results with lattice data [T7] and 



The T = onset of diquark condensation at jU c (0) = 
m T /2 as an exact result is built-in in xPT and the 0(6) 
linear sigma model as discussed in Sec. II B Therefore, 
it also holds for the screening mass of the pion from the 
bosonic potential in our adapted no-sea approximation, 
which reduces to the linear sigma model at T = by 
definition. They both go beyond the leading order chiral 
perturbation theory in that they include effects of a finite 
sigma meson mass. The linear sigma model expressions 
for the T = condensates are [Ml 




for fi < [i c 
for jj, > ji c 



1 - 



25 



for n < \i c 
for ji > /i c 



(25) 



where x = 2///m,r and y = m a /m 7T . The only difference 
between these and the x?T result [7J is the y-dependent 
term in the diquark condensate which reduces to the xPT 
formula for y — > 00. Note that the chiral condensate does 
not depend on the sigma meson mass whereas the diquark 
condensate does via y which explains the variations of the 
diquark condensates at large /i in Fig. [3j 

Beyond the no-sea approximation, one needs to distin- 
guish between screening and pole mass. Only the latter 
agrees with the onset baryon chemical potential at the 
mean-field level, in general. We will discuss this in more 
detail in the next subsection. 

The overall agreement of our T — results with the 
existing literature in Fig. [3] is very reassuring. In par- 
ticular, there is no dependence on the chemical potential 
for (i < fi c in accordance with the Silver Blaze property 
as also discussed more in Sec. MI Cl below. 

Common to all studies, the chiral condensate decreases 
with increasing in the diquark condensation phase 
above /x c . As a result, the quark- meson- model-like 1 st 
order transition line and CEP at /1 around 2.5to t are 
completely gone, as seen also in Fig. [I] 

Fig. [4] shows the phase diagram of two QMD model 
mean-field calculations in comparison with the result of 
an NJL model recalculation following Ref . [5S] . For small 
chemical potentials and small temperatures one finds the 
chirally broken phase with vanishing diquark condensate 
crossing over with increasing temperature to the phase 
with asymptotically restored chiral symmetry as usual. 
In addition, there is the diquark condensation phase for 
fi > /i c (T) characterized by non-zero baryon density. 

The NJL model calculations show a continuous diquark 
condensation transition throughout the whole phase di- 
agram [28 . In contrast, in our mean-field QMD model 
results here, the second order line /x c (T) ends in a tri- 
critical point from where on towards larger chemical po- 
tentials it becomes a 1 st order transition. Such a behav- 
ior was also predicted at next-to-leading order xPT [9], 
where it was concluded that this tricritical point occurred 
at n c « 0.57?^ and T w 220 MeV. The temperature 




NJL 

QMD MF A - 600 MeV 
QMD MF A - MeV 



IJ [««] 



FIG. 4. Mean field phase diagrams: QMD model (with and 
without vacuum term) vs. NJL model (with parameters from 
[26]); dashed: chiral crossover; solid: second order transition; 
thick solid black: first order transition; diamonds: tricritical 
points as predicted in [9]. 



C. Pole versus screening masses 

The definition of meson/diquark masses is absolutely 
crucial already at mean-field level to get results which 
are consistent with the Silver Blaze property. As soon 
as one goes beyond the no sea approximation there is an 
important distinction between the screening mass at zero 
external momentum, which is just determined by the cur- 
vature of the effective potential, and the pole mass which 
takes into account a non-vanishing external momentum. 
The latter is the natural choice in NJL model calculations 
but has so far not been taken into account in QM model 
studies. Both definitions coincide for massless particles 
and in the no sea approximation. 

The meson/diquark pole masses are defined as the ze- 
roes of the determinant of their inverse propagator, 



ri? ) (p) = r t w( P ) ii + n tf ( P ) 



,(2 



(26) 



where, for the pions, 

.(2) 



r;r(p)« = (P 2 - ™ 2 + M 2 )5 l3 , i,j = 2,3,4, (27) 



range is quite comparable here, but its precise location 
depends on the value for the vacuum-term cutoff A and 
moves towards larger chemical potentials with increasing 
A. Moreover, as we will see below, the first-order transi- 
tion will be washed out by fluctuations which will make it 
second-order again, and it is hence a mean- field artifact. 



with 

A3 



m 

.2 



= Xv 2 from the tree-level potential (10 1 and 
= er"+G? 2 . When the diquark condensate d is non-zero, 
we choose it to lie in the 05 = Re A direction without 
loss, the sigma meson mixes with the scalar diquark pair 
already at tree-level through the 0(6) linear sigma model 
potential (10). In the 3-dimensional subspace of sigma 
and diquarks in the real basis 05 = Re A, 06 = Im A, 
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+ A0 2 + 2Xa 2 
2\ad 




p 2 - An 2 



2Xad 
m 2 + A0 2 



2Ad 2 




(28) 



The RPA polarization functions are obtained from evalu- 
ating the fermion-loop integrals with external momentum 
p in the usual way, 



IIy(p) = Tr, 



dT 



(2) 



0MF 



Gmf(p - 



9r (2) 



<f>MF 



Gmf(q) 



(29) 

where Gmf = (rpP|^ MF ^ . They can be found in the 

NJL model literature 28 , 39 , 40 , originally from the two- 
flavor three-color standard NJL model with isospin chem- 
ical potential and pion condensation. Since the available 
expressions are either incomplete or at variance with our 
computations, we have recomputed them and list the 
complete explicit expressions for these polarization func- 
tions as a convenience to the reader in Appendix [C| 

To find the pole masses we use p = (— iw,0) and the 
somewhat sloppy notations r' 2 )(u;) = T^(p = (— iw,0)), 



n(w) = = (-iw, 0)), to solve 



detr (2) (w) = , for w = ra fc , k = !,...&. (30) 



The polarization integrals are ultraviolet divergent. As 
before, we use a spatial momentum cutoff A for the 
temperature-independent contributions. Making the T- 
dependence explicit, we may thus write, 



n rcs (w,T) = U th {cj,T) +nX ac (w) , where 

n th (w, t) = n(w, t) - n(w, o) 



(31) 



is ultraviolet finite. In the normal phase with d — 0, i.e., 
for fi below the onset of diquark condensation at /i c (T), 
the polarization integrals are diagonal in the basis where 
A = 4>5 + i(f>6 an d A* = 05 — i(f>6. The polarization matrix 



9 



n(w, T) from App. [C]is then diagonal with entries 



n CT (w,T ) =16iV c5 /t^3 2 , 2 — — 



n±(a;,T) =4iV c5 : 



d 3 q e q {l-N q (e-)-N q (e+)) 



uj 2 - 4e 2 



(2tt) 2 



2ef uj + 2ei 



(32) 



with = y^T^Vi n, n± for A, A* and Polyakov 
loop enhanced quark/antiquark occupation numbers 



N q (E) = N q (E;T,$) 



1 + $er 



1 + 2$er 



2E 3 



(33) 



which simplify to the Fermi-Dirac distribution for $ = 1. 

As usual, these expressions are obtained from analyti- 
cally continuing the results for imaginary discrete values 
uj = \2uTn corresponding to the Matsubara frequencies 
in imaginary time. To make the continuation unique, one 
usually assumes in addition that the polarization func- 
tions are well behaved at complex infinity with cuts only 
along the real axis. Then, it follows that the correspond- 
ingly continued expressions for finite spatial momenta, 
H(p = (—iu>,p)) are non-analytic at the origin in mo- 
mentum space, with different limits for ui — > at p = 
or for \p\ — > at uj — 0. The first limit yields a plasmon 
mass, associated with the damping of plasma oscillations, 
while the second is the one that yields the correct finite 
temperature screening masses [54) . Here, in the normal 
phase the two differ only for the sigma meson, by the 
n = contribution proportional to 5 U fi in the equation 
for n <T (o;,T), which can be obtained from the expression 
for H a {p = (— \u),p)) in Ref. [28] with the additional pre- 
scription to set uj = first and then take \p\ — >• 0. 

The corresponding extra contributions for w = in 
the diquark condensation phase are also given in App.[C| 
None of them are really needed here. In particular, the 
^,oII (T) contributions vanish for T — > 0, but they as- 
sure that the screening masses extracted from the propa- 
gators agree with the corresponding ones from the effec- 
tive potential also at finite temperature, see below. 

Setting $ = 1 and dismissing the temperature depen- 
dent contributions S^filL (T), the polarization functions 
agree with the ones from Ref. [351 HO] for baryon chemical 
potential /ib — and isospin chemical potential fii = 2/i, 
where II ^ and II -t belong to neutral and charged pions, 
respectively. 

The RPA pole masses in the quark-meson-diquark 
model in the normal phase are then simply given by the 



solutions of the following equations, 

uj 2 = -m 2 + 3\a 2 +II CT (w,T) 
uj 2 = -m 2 + Act 2 +U 7r (uj,T) (34) 
(w±2/i) 2 = -m 2 + Act 2 +n±(w,T) 



rn 



m± 



If we use the mean-field expression in Eq. ( 19 1 for the 



(2tt) 3 

d 3 q (l-2N q (eJ) l-2JV,(e±)\ Q q (T, = - 



fermionic pressure with chiral and diquark condensates 
in the form, 



4T /(|^I>( 2 ( C0Sh 



T 



one immediately verifies that the polarization functions 
for external momentum p = 0, corresponding to the limit 
\p\ — > at uj = in the imaginary time formalism, 



n CT (o,T) 

IMO.T) = 2 



-4cr 2 



d 2 



da 2 



d=0 



d=0 

d 



2 \2 



n q (T,n) 



d=0 



(36) 



n+(0,T) = n_(0,T) = 2 — Sl q {T,n) 



d=Q 



This shows explicitly that the screening masses, defined 
by these derivatives of the effective potential, are ob- 



tained as the constant contributions in Eqs. (34) for 
0, 



UJ 



m s £ 2 — —m 2 



\a 2 + U a (0,T), 

m sc2 _ _ m 2 + A(T 2 +Un (0,T), 

m s l 2 = -4fi 2 - ni 2 + Act 2 + U ± (0, T) 



(37) 



This is true at all temperatures in the normal phase. 
Note also that because n+(0, T) = II_(0, T), the baryon 
chemical potential fiB = 2/i never splits the diquark and 
antidiquark screening masses, m s ^(T,fj,) = m^(T, /z). 

At any temperature we furthermore verify for // = 
that H n (uj,T) = II±(w,T), i.e., pion and diquark masses 
are degenerate as they must from 5*0(5) symmetry. 
Moreover, the gap equation for the chiral condensate re- 
duces in the chiral limit c — > to the condition for mass- 
less pions, as usual, and both these observations hold for 
screening and pole masses, likewise. 

Finally, but maybe most importantly, the gap equation 
for the diquark condensate reads 



d_ 
~d~d" 



A<7 2 -4/i 2 



0, (38) 



and the critical line /i c (T) is defined by the condition that 
the terms in brackets vanish for d = so that a second 
zero develops there. This is equivalent to the diquark 
pole masses being m_ = and m + = 4/i. While their 
screening masses m± both vanish at fi = fi c , for the pion 
and diquark pole masses we have the general exact zero- 
temperature relation 



n±(w,o) = iUw±2 M ,o) 



m± = ni-n ± 2/i , (39) 
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in the normal phase where = m,^ remains indepen- 
dent of /i until 2/i — m Wj o &s required by the Silver Blaze 
property. 

In contrast, the same relation entails for the degenerate 
diquark screening masses (II W is an even function of w), 



-,sc2 



^sc2 



V + n ff (2 / i,o)-n, r (o,o), (40) 



which reiterates that diquark and pion screening masses 
are also degenerate at fi = 0, but that both diquark 
screening masses m± vanish as 2/i approaches the (/x- 
independent) pion pole mass m, from below. 

The bottom line is that the onset of diquark conden- 
sation at hb = 2/i = 2/z c (0), whatever the screening 
mass may be, defines the physical zero-temperature pion 
mass. We will make use of this property to fix the pion 
mass in the RG calculation, where the calculation of the 
pole mass is more involved. 

In the diquark-condensation phase the sigma meson 
mixes with the two diquark modes, i.e., the respective 
masses have to be determined from the zeroes of the 
determinant of the corresponding 3x3 submatrix in 
r( 2 )(w). As in the NJL model [Ml ES] one can verify 
further exact results from the mass formulas at T = 0. 
Also in the QMD model at mean-field/RPA level the 
in-medium pion pole-mass is equal to m w = 2/i above 
the onset of diquark condensation at 2fi — m W fi. More- 
over, one verifies explicitly that one of the three modes 
in the diquark/sigma sector remains exactly massless in 
the superfuid phase, also at finite temperature. This is of 
course the Goldstone boson corresponding to the spon- 
taneously broken U(1)b baryon number. Another one 
becomes degenerate with the pions for large values of the 
chemical potential, eventually, reflecting the restoration 
of chiral symmetry. They combine into an SO (4) multi- 
plet as the chiral condensate vanishes for large chemical 
potentials. 




FIG. 5. Pole-mass spectrum at T = 0: mean-field/RPA QMD 
model results (with vacuum-term cutoff A = 600 MeV), for 
comparison also shown without the effect of diquark/sigma- 
meson mixing in the superfluid phase, vs. linear sigma model. 



This is all nicely reflected in the numerical results 
shown in Fig. [5j As the RPA pole-mass formulas im- 
ply, the meson masses stay constant in the normal phase 
whereas the diquark masses are split up from the con- 
stant m,B — by the terms ±/ib due to their coupling 
to the baryon chemical potential \ib = 2/x. 

The diquark and sigma masses in the phase of diquark 
condensation show a considerable dependence on the in- 
clusion of the vacuum term. This can be seen, for exam- 
ple, by comparing the QMD model results with vacuum- 
term cutoff A = 600 MeV to those from the linear sigma 
model, which are identical to the ones in the no-sea ap- 
proximation (A = 0). In the linear sigma model, the 
pole masses can simply be calculated from the curvature 
of the potential. In the normal phase they are simply 
given by the expected constant m x = TOtt.Oj m a = "V.cb 
and m± = m,j ± 2/i. In the phase with diquark con- 
densation (/i > /i c ), on the other hand, we obtain for the 
linear-sigma model masses, 



m, =2/i m^j = 



m a/A2 =^ (-3 + y 2 + 28x 2 ± (3(1 - y 2 )/ 



2\ / 2 
1 'X 



+ (-3 + y 2 ) 2 + 40(-3 + y 2 )x 2 + 400x 4 ) 3 



(41) 



with x = 2/1/771^ and y — m^/m^ as in Eq. (25). 

In addition, Fig. [5] also shows results of a mass cal- 
culation where the mixing terms in the sigma/diquark 
sector were neglected. In those results we have explic- 
itly set the off-diagonal tree-level mixing term 2Xad in 
Eq. (28) to zero, and used IT^a = TlaA* — for the po- 
larization functions of App.[C]in the superfluid phase. As 
illustrated in Fig. [5] while the pole masses in the no-sea- 
approximation are very close to the xPT result, the only 
difference is due to the finite sigma mass in Eq. (41), 



the more realistic ones with sufficiently large vacuum- 
term cutoff are closer to those without any mixing in the 
crossover region at intermediate chemical potentials. 



IV. FRG FLOW EQUATIONS 

Quantum and thermal fluctuations arc of utmost im- 
portance in particular near phase transitions and are con- 
veniently included within the framework of the functional 
renormalization group (FRG) [551461] . In this work we 
employ a Wilsonian RG version and investigate the flow 
equation for the effective average action pioneered by 
Wetterich [55] . The central object in this approach is 
the renormalization scale k dependent effective average 
action [$] , where $ generically represents the set of all 
quantum fields of the theory. The effective average action 
interpolates between the microscopic classical action at 
some ultraviolet (UV) cutoff scale k = A, at which fluctu- 
ations of essentially all momentum modes are suppressed, 
and the effective action of the full quantum theory in the 



11 



infrared (IR), for k — > 0, which then includes all quan- 
tum and thermal fluctuations. The scale-dependence is 
described by the Wetterich flow equation, 

d t T k = kd k T k [^>] = ^Tr { d t R k (T^ + Rk)- 1 } , (42) 

which involves a momentum- and scale-dependent reg- 
ulator R k , whose precise form is not fixed but leaves a 
considerable flexibility. The role of the regulator R k is 
to suppress the fluctuations of modes with momenta be- 
low the renormalization scale k, and the flow equation is 

(2) 

UV as well as IR finite. T k [$] are the second functional 
derivatives of the effective average action with respect 
to all the fields at scale k. The functional trace repre- 
sents a one-loop integration typically evaluated in mo- 
mentum space and includes the sum over all fields and 
their internal and space-time indices as well, with stan- 
dard modifications for fermionic fields. It contains the 
full field and fc-dependent propagators of the regulated 

(2) 

theory with cutoff R k , the inverse of TV [$] + R k . In 
order to solve the flow equation an initial microscopic 
action S — T k= A at some UV scale A has to be speci- 
fied. Truncating the effective action to a specific form, 
the functional equation can be converted into a closed 
set of (integro-)differential equations, but will in general 
also introduce some regulator dependence in the flow. 
The choice of an optimized regulator minimizes this reg- 
ulator dependence for physical observables. As bosonic 
(fermionic) regulators R kt B [Rk,F) we choose 



( fk* \ (43) 
R k . F (p) = -ip- f - 1 j 9(k 2 - f), 

which are three-momentum analogues of the optimized 
Litim regulators |63j . With this choice the three- 



momentum integration becomes trivial and the remain- 
ing Matsubara sums can be evaluated analytically. Fur- 
thermore, this choice leaves the semilocal [/(l)-symmctry 
of the Lagrangian unaffected, analogous to pH], where 
the chemical potential acts like the zero-component of an 
Abelian gauge field. In addition, these regulators have 
their precise counterparts in specific three-momentum 
regulators for proper-time flows which then lead to iden- 
tical flow equations, c.f. Appendix|E] On the other hand, 
especially in a relativistic system the fact that the zero- 
component of the momentum is not regulated can poten- 
tially be problematic. 

Our ansatz for the effective average action in leading- 
order derivative expansion, where all wave-function 
renormalization factors are neglected and only the scale- 
dependent effective potential U k is taken into account, 
reads 



r fc = 



d 4 x£pQMD 



(44) 



V+ca^U k 



This means that we use £pqmd from Eq. ( 14 ) , but re- 



place the meson/diquark potential V{4>) of the 0(6) lin- 
ear sigma model from Eq. ( 10 1 therein by U k — ca. The 
explicit symmetry breaking term —ca does not affect the 
flow and is thus not part of U k but added after the RG 
evolution to the full effective potential again. At p = 0, 
the scale-dependent U k then only depends on the modu- 
lus of 4> — (a, if, Re A, Im A) T . At non-vanishing chemi- 
cal potential, however, we only have 5*0(4) x 5*0(2) sym- 
metry and must therefore allow it to depend on two in- 
variants, i.e., U k = U k (p 2 ,d 2 ) where p 2 = a 2 + if 2 , and 
d 2 = |A| 2 as before. For p — > we recover 50(6) in- 
variance, of course, so that U k then depends only on the 
combination <fi 2 = p 2 + d 2 again. 

With the constant field configurations cr = p, n = 0, 
Re A = d, ImA = we obtain for the bosonic second 
functional derivative of the effective action 



(p 1 



r (2) _ 
1 k,B — 



h 2U k , p 










- 2U k , 







p 






+ 2t4„ 






p 






- 2U k , p + 4p 2 U k , 
4pdU k , pd 








p 2 + 2U k4 + Ad 2 U k4d 



Ap? 








-4/xp 
■ 2U k , d - 



(45) 



where we have introduced short-hand notations for the 
derivatives of the potential with respect to the fields de- 
fined as U k ^ = dU k /dd 2 , U k p = dU k /dp 2 and later we 
will also use U k ^ = dU k /d(f> 2 . Higher order derivatives 
are labeled with higher order indices accordingly, e.g., 
U k , p d = d 2 U k /dp 2 dd 2 . The alert reader will have no- 



ticed that Eq. (45) agrees with Eqs. (27 1 and (281 upon 



working out these derivatives, if we replace U k back to 



X((p 2 — v 2 ) 2 /4 which is what we use at k = A. 
In the fermionic sector we find analogously 



.(2) 
k,F 



-i|»-iao7 +SP-7 M 



-97 d 



91 d 

-i^-iao7° '+gp+7 n M 



^2x2 



(46) 

To these expressions we add the respective regulators 
in Eq. ( 43 1 before they are being inverted and inserted 



into the Wetterich equation Eq. (42), in order to ob- 
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tain the flow equation for the effective potential. Re- 
placing the zero-components po of the momenta by peri- 
odic (antiperiodic) Matsubara frequencies ui n = 2irnT 
(y n = (2n + l)7rT), and upon performing the spatial 



momentum integrations, the corresponding bosonic and 
fermionic contributions to the flow for the effective po- 
tential are then given by the following Matsubara sums, 



dtUh.B = 



k 5 T 



V 



d t U K 



k - B 6ir 2 ^ 

nez 

8k 5 T 



a 2 (u 2 l ) 2 + aiw, 2 , + a Q 



Ljl + k 2 +2U k>p K)3 +/ 3 2 ( w 2 )2+/3lW 2 +A) 

(u n +a ) 2 + k 2 + g 2 <j) 2 - t i 2 



^ ^ z {{v n + a Y + Ef)((u n + a,Y + Ef) ' 



(47) 
(48) 



where we have introduced E^ — -J g 2 d 2 + (e k ± /i) 2 , and inverse of the 3x3 submatrix of the bosonic 2-point func- 
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£fc = \/k 2 + g 2 p 2 analogous to the notations of Sec 
The numerator of the second term in the bosonic flow, 
Eq. (47), is a quadratic polynomial in u 2 with three co- 



efficient functions ai, while the denominator is a cubic 
polynomial in standard form, with coefficient functions 
Pi and leading coefficient f3 3 — 1. These coefficient func- 
tions ai and Pi depend on renormalization scale, chemi- 
cal potential, fields and the derivatives of the potential. 
They are obtained straightforwardly from the trace of the 



tion in Eq. ( 45 ) corresponding to the sigma and diquark 



directions in field space, and they are listed explicitly in 
Appendix ID] for completeness. With the roots of the de- 



nominator which we denote as i 



'n,0 



1, 3, we 



can evaluate all Matsubara sums analytically by virtue 
of the residue theorem in a standard way. 

Hence, the final flow equation for the effective poten- 
tial of the PQMD model is the sum of the bosonic and 
fermionic flow and reads explicitly, 



d t U k = 



12tt 2 







— coth | 




K 2T) 



E a 2 zf - aiz 2 + a Q 1 / Zj \ 

, , {z 2 +1 -z 2 ){z 2 +2 -z 2 )zi\2T) 



1 ± 



yjk 2 + g 2 p 2 



(l-2JV g (^;T ) $))| , 



(49) 



where E£ 



sjk 2 + 2U k , P , and N q (E;T,<f>) are the 
Polyakov loop enhanced quark occupation numbers from 
Eq. (33). Without a diquark condensate, i.e. by setting 



explicitly A = 0, we can write down an S'0(6)-symmetric 
flow equation for U k ((j)), if we set Vh,4> = U k , P — U k ,d- 
Eq. ( 49 1 then reduces to the more familiar looking form, 



dtU k = 



I2ir 2 



El 



coth 
1G 



El 
2T 



1 

EZ 



coth 



E% 
2T 



El 



coth 



El - 2m 
2T 



1 - N q (e fc - M5 T, $) - N q {e k + /x; T, $) 



1 , (El 
— coth ' K 

El 



2n 



2T 



(50) 



with single-particle energies for mesons/diquarks El = 
yjk 2 + 2f7 fc ,0 and sigma E% = y/k 2 + 2E/ M + ±<t> 2 U ky<t><t> . 
Except for the change in the number of active degrees 
of freedom contributing to this flow, and the isospin-like 
chemical potential coupling to one pseudo-Goldstone bo- 
son pair, the SO(6) symmetric flow equation here is en- 
tirely analogous the one of the PQM model in the three- 
color case, see e.g., [48, 64, 65 . For the three-color PQM 



model with isospin chemical potential one must allow for 
pion condensation, however, and then arrives at a flow 
equation [41 analogous to our Eq. (49). 



In the following sections we present numerical solutions 
to the flow equation (49). The full effective potential 



depends in general on three condensates which hampers 
its numerical solution enormously. In order to proceed 
we restrict ourselves in this work to the two-color QMD 
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FIG. 6. Fit to chiral condensate for critical exponent /3. 



FIG. 7. Fit to chiral condensate for critical exponent S. 



model and neglect the influence of the Polyakov-loop by 
setting $ = 1 in the flow and postpone the full PQMD 
model solution for a later analysis. For the first time we 
generalize the one-dimensional grid solution technique to 
two dimensions. Details of the numerical procedure and 
the parameter fixing can be found in Appendix [X) 



A. Critical exponents f3 and 5 

Without diquark condensation for vanishing chemical 
potential, and Polyakov-loop variable $ = 1, the 5*0(6) 
symmetric flow in Eq. ( 50 ) further simplifies to 



rameter c, 



d t U k 



12tt 2 



El 



coth 



El 
2T 



1 

El 



■ coth 



El 
2T 



16 



tanh ( ) 
e k \2T) 



(51) 



At n = the diquarks are degenerate with the pions 
which leaves us with the N c — 2 analogue of the familiar 
three-color QM model flow equation [SSI EZ] except that 
there are now five pseudo-Goldstone bosons instead of 
the usual three pions. 

The study of 0(4) universality and scaling in the three- 
color QM model has a long history by now [58HT2"] . Here 
we can analogously check the symmetry breaking pat- 
terns discussed in Section |II A| by computing the cor- 
responding critical exponents. As discussed there, for 
/i = rn — 0, the SU(4) ~ 5*0(6) dynamically breaks 
down to Sp(2) ~ 50(5) so that we expect a finite tem- 
perature phase transition in the three-dimensional 0(6) 
universality class. The critical exponent j3 can be ex- 
tracted from the dependence of the chiral condensate on 
the reduced temperature t = (T — T c )/T c in the chiral 
limit, whereas the exponent 5 governs the dependence 
of the chiral condensate at T c on the quark mass m q or 
correspondingly on the explicit symmetry-breaking pa- 



(qq) T ~ (-tf, (qq), 



(52) 



With the usual two-exponent scaling all other critical 
exponents are then obtained from these two. Here we 
find critical exponents /3 = 0.4318(4) and S = 5.08(8) 



from the solution of the lei flow equation (51 1 via the 
Taylor expansion method. The given errors are statis- 
tical errors extracted from the fit. The corresponding 
fits are shown in Figs. [6] and [7] respectively. Literature 
values for these exponents obtained from Monte-Carlo 
simulations are given by /? = 0.425(2) and 5 = 4.77(2) 
[73] . At the leading order derivative expansion employed 
here, we should not expect to reproduce these values, 
however. The more appropriate benchmark here should 
be the functional renormalization group result for the 
0(6) model in leading order derivative expansion [74j . 
In absence of wave-function renormalizations there is no 
anomalous dimension for the fields and their critical ex- 
ponent therefore vanishes, 77 = 0. Then the hyperscaling 
relations, 

immediately entail that S = 5 and /3 = v/2 in d = 3 
dimensions. With the correlation-length critical expo- 
nent v = 0.863076 from Ref. [71] this corresponds to 
f3 = 0.4315, and both our values are in agreement with 
these two within our errors. 



B. Phase diagram without diquark fluctuations 

Before we discuss the solutions to the full flow equa- 
tion (49) for the effective potential with fluctuations of 
included, 



both condensates included, it might be instructive to 
illustrate the influence of fluctuations on the standard 
quark-meson-model-like phase diagram without baryonic 
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degrees of freedom. The phase diagram obtained from 
our solutions to the S'0(6)-symmetric flow equation (50) 



in the (T, /i)-plane is compared to the mean-field results 
from Sec. |III A| in Fig. [8 The mean-field solutions there 
were obtained from Eq. (19j with d 2 = |A| 2 = and with 
a vacuum term cut off at A = 600 MeV which is suffi- 
ciently large for a reasonable comparison, see Sec. |III A| 
The parameters were chosen so as to match the fi = chi- 
ral transition temperatures (rather than the sigma mass) 
in addition to pion decay constant and pion mass, as ex- 
plained in Appendix |A") 

Again, the resulting phase diagram with fluctuations 
shows the typical form of the QM model for N c = 3. 
It has a critical endpoint at /i « 2.5m 7r as compared to 
that at fx « 2.8/71^ in the mean-field calculation. The 
dotted chiral-crossover lines are again simply the half- 
value curves of the chiral condensate. Except for the shift 
of the critical endpoint and the crossover-line in the range 
of chemical potentials between 2.5 and 2.8 m T , the chiral- 
crossover line from the RG solution closely follows the 



mean-field result. As discussed in Sec. Ill B already, they 
are both equally wrong for baryon chemical potentials 
of the order of the baryon mass and above, as we have 
neglected the essential dynamics of the baryonic diquarks 
in the superfluid phase. 
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FIG. 8. QMD model phase diagram without baryonic degrees 
of freedom (d = 0): Id RG result vs. MF (A = 600 MeV). 

The advantage here as compared to the three-color 
case is that it is much more straightforward to include 
baryonic degrees of freedom, by properly including the 
diquarks, both at mean-field level and in the full flow 
equation ( 49 1 . All we have to do is to solve this equa- 



tion in the higher-dimensional space of invariants in the 
fields representing both possible order parameters and 



their fluctuations at the same time. 



C. FRG pole mass and flow of the 2-point function 

The fact that the Silver Blaze property links the on- 
set of diquark condensation to the zero-temperature pion 
mass represents a strong constraint which has to hold in 
the (P)QMD model with fluctuations also. As we have 
already illustrated at mean-field level in Section |III C[ a 
proper definition of meson and baryon masses is abso- 
lutely crucial for this exact property of the theory. In 
particular, we have seen explicitly that the behavior of 
widely used screening masses is unphysical in this regard. 
In the mean-field calculations, the difference between the 
pion's pole and screening masses at zero temperature and 
with identical model parameters can be as much as 30%, 
for example. Adjusting the parameters to the more phys- 
ical pole mass instead of the common procedure in these 
models has a considerable influence on the results. So 
this is more than an academic exercise. 

Therefore, we propose a simple procedure to obtain 
pole masses suitable for the FRG framework: As an ex- 
tension to the flow equation for the effective potential 
in the leading order derivative expansion, we solve the 
flow equations for the 2-point functions of mesons and 
diquarks using the field and scale dependent but momen- 
tum independent 3- and 4-point vertices obtained from 
the effective potential. This ensures that the flow equa- 
tions for the 2-point functions at zero momentum reduce 
to those for the mass terms in the effective potential and 
that the screening masses obtained from the flows of 2- 
point functions and effective potential are therefore the 
same by construction. This truncation for the flow of the 
2-point functions most naturally extends that of the ef- 
fective potential, and thus provides a simple alternative 
to other approaches such as the BMW approximation 
[75] where the same relation with the effective potential 
typically arises as an additional requirement. 

In this section we outline the derivation and solution 
methods for the flow equation of the pion 2-point function 

(2) 

rW at T = /i = as an example which will allows us 
to define a pion pole mass in the vacuum. We consider 
Nf flavors of quark with N c colors coupled to an O(N)- 
symmctric bosonic sector for combinations of Nf , N c and 
N where this is possible. 

The flow equation for the (field dependent) 2-point 
function is given by the second functional derivative of 
the original flow equation (42) which in our case is, 



d t rg< 2 W) = J ^a f i? B ( g ) H (G B ( g ) im r^G B (g- P ) r , r r^ 3) G B ( g ) sfc - \G B { q ) lm vt%G B { q ) nk ) ^ 

-2trd t R F {q)(G F { q )Tf^G F { q - P )Yf^G F { q ) - \g f { q )Tf^ G F {q)) , 
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where G B (p) 

(p;0)- 



.(2,0) 



(T^(p;<p)+RB(p)) 1 and G F (p)_ 



■R%(p)) ■ To solve the flow equation (54) 



one needs 3 rd and 4 th derivatives of I\. [<fi] which we de- 
note generically by rj c ™' m ' 1 where the first(second) super- 
script counts the number of fermionic(bosonic) deriva- 
tives. In order for the limit p — > to be consistent with 
the truncation used for the flow equation of the effective 
potential we obtain those higher n-point vertex functions 
from the same scale-dependent effective action. In lead- 
ing order derivative expansion, Eqs. (44), (49 1 or (51 1 for 
fi = 0, the 3- and 4-point functions are then momentum 
independent, and the only dependence on the external 
momentum comes from the propagators themselves. For 
convenience we choose coordinates <f>i = <f>6n, i.e., a = <f> 
and the others zero, so that one has explicitly for the 
quark-boson vertices with constant Yukawa couplings, 



r 



(2,1) _ 



ri(2,l) . o „ 



(2,2) = 
ij 







(55) 



The three and four-boson vertices are extracted from the 
respective derivatives of the fc-dependent effective poten- 



tial Uk ■ Here this is a function of (j> 2 and we simply write 
U'f!((f> 2 ) = Uk,4><t> etc. instead of our index notation above, 



,(0,3) 



= # U'l (SijSmi + Sj m 8ii + S im Sji) 



(56) 



,(0,4) 



— 4 Ufa ($ij$mn ^in^jm &jn$im) 



+ Sj n 5n5 m i + SinSjiSmi) 



(57) 



For the calculation of the boson masses we use their rest 
frame, setting the spatial external momentum p = 
in Eq. (54). In this frame the spatial momentum in- 



tegrals with the optimized regulators are still trivially 
performed. Evaluating the flow equation (54) one ob- 
tains after analytically continuing pq = — iu;, and with 
the sa me not ations r (0 > 2) (w) = T (0 ' 2) (p = (-iw,6);</>) as 
in Sec. |IIIC| above, 



ftr£2( w ) 



(N + \)U'i , 2UZ{Ef-Ef){{El+EZf{Ef + E° k EZ + Ef)-{Ef 



Er 



U£ + 2fU, 
El 



(3) 



8N f N c g 2 (4el 




lo 2 Y 



(58) 



We set T 



(0,2) 



-(J 



1 ) = -uj 2 + X{<j ) 2 -v 2 ) 



at the UV scale k = A, and obtain the pion pole mass 
«T"7r,poic for k — > from the condition 



.(0,2) 

k— 0,7T7r 



(«4, P o 



0. 



(59) 



evaluated at the minimum of the full effective potential. 
For vanishing external momentum the two-point function 
can equally be obtained from the second derivative of 
the effective potential. Indeed, one verifies that the flow 



equation ( 58 ) obeys the consistency condition 



a r (0,2) 



(0) 



cJ 2 



N- 1 d^d-ir 



-d t Uk 



2-^dtU k . 



(60) 



This implies that if we calculate T^f^^ by integrating 



the flow equation (58) for cj — 0, the mass defined as 



< 2 = r^(o) 



(0,2) 



(61) 



correspondingly, simply represents the same screening 
mass as obtained from the curvature of the effective po- 
tential at its minimum, which is usually considered in 
QM model calculations within the FRG framework. 
The flow equation (58) can be solved via a Taylor 



expansion method around a scale dependent expansion 



point for both the effective potential and the two-point 
function, or on a grid in field space. In order to maintain 
the relation in Eq. ( 60 1 also in the numerical calculations 
based on Taylor expansions in (f> 2 , one has to use one ex- 
pansion order less for the 2-point function than for the 
effective potential. In this way we can compute an esti- 
mate of the pion pole mass from a given UV potential. 

Table [T] shows a comparison of screening and pole 
masses as obtained from the Taylor and grid methods. 
All calculations here were performed at T = fi = 0. As 
explained in Sec. |III C| and Appendix[A] we have adjusted 
the start parameters for the flow in our two-dimensional 
grid code to fix the onset of diquark condensation to oc- 
cur at 2/j, c w 138 MeV which defines the physical pion 
mass in the normal phase. The exact same parameters 
were used to obtain the UV forms of effective potential 
and inverse propagators for the one-dimensional Taylor 
expansion method at y, = 0. The results from one- and 
two-dimensional grid computations at [i = are indistin- 
guishable at this level of accuracy, as are the screening 
masses from Eq. (61) and from the effective potential. 
The slight deviations in / x and the masses in Tabic [I] be- 
tween the grid and Taylor methods are an indication of 
the small residual uncertainties. 

With the onset at half the physical pion mass fixed, 
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FIG. 9. Zero temperature condensates from full flow com- 
pared to mean- field results (and the lattice data from [17]). 



we then observe that the standard screening masses gen- 
erally overestimate the pion mass by about 30%. In 
contrast, our pion pole mass estimates based on solving 
Eqs. (58 1 and (59 1 lie within 11%, but they are smaller 



than the physical one. 

The extrapolation from zero pion momentum in the 
leading order derivative expansion to the pion pole in 
the propagator from our consistent truncation scheme 
appears to be too large, so that it overcompensates the 
difference between onset and screening mass. In order 
to disentangle the effect of bosonic and fermionic con- 
tributions to the flow equation ( 58 1 for the pion 2-point 



function, we have also solved this equation with u) = in 
the bosonic and in the fermionic parts, separately. The 
resulting pole masses are denoted by m^poio, f orm . on i y 



method 


quantity 


value [Me V] 


Grid 


u 




76.0 








178.8 








551.7 




2/i c 




137.8 




^"7T,polc 




122.45 




T~H>tt .pole, ferm. 


only 


124.9 




771-77, pole, bos. 


only 


171.6 


Taylor 


/* 




75.0 




TYl-Tr ,scr 




180.0 








550.8 




^"7r,polc 




122.6 




^^71", pole, ferm. 


only 


125.0 




f^7r,pole, bos. 


only 


172.6 



TABLE I. Comparison of RG screening vs. pole masses; 
'ferm. only' ('bos. only') refers to maintaining only the con- 
stant lj = contributions in the bosonic (fermionic) contri- 
bution to the flow of the pion 2-point function, Eq. A58\. 



and TOTr.poic, bos. only hi Table [TJ respectively. Both con- 
tributions reduce the screening masses, but the fermions 
clearly generate the dominant effect. This suggests that 
one might have to go beyond the leading-order derivative 
expansion employed here and allow for an RG flow of the 
Yukawa couplings by including field renormalizations and 
anomalous dimensions [76j . 



D. Phase diagram of the QMD model for two-color 
QCD with mesonic and baryonic fluctuations 

In Fig. [9] we show once more the dependence of the 
chiral and diquark condensates on the chemical potential 
at zero temperature as in Fig. [3j but this time with in- 
cluding our results from the full RG solution to Eq. ( 49 ) 
on a two-dimensional grid in field space with $ = 1. 

The final effect of baryonic-diquark degrees of freedom 
is illustrated in Fig. [10] where we compare the phase di- 
agram from the one-dimensional RG flow solution to the 
50 (6)-symmetric equation 



(|50J) from Sec. |IVB| and that 
from the full two-dimensional one for an effective poten- 
tial with the reduced 5*0(4) x 50(2) symmetry. 

This clearly illustrates the effect of the competing dy- 
namics of the collective mesonic and baryonic fluctua- 
tions. As before, the dashed lines in Fig. [lO] indicate 
the chiral crossover by tracing the half- value of the chiral 
condensate. Both, the one and the two-dimensional re- 
sults agree for quark-chemical potentials near zero. The 
crossover in this region leads to mesonic freeze-out as 
usual, and the results are unambigously determined by 



the 0(6) symmetry breaking pattern, see Sec. IV A Al- 
lowing additional interactions with lower symmetry has 
no effect on the flow here. 

Once the quark-chemical potential approaches half the 
baryon mass, corresponding to ms/7V c , however, the 
rapidly increasing baryon density equally rapidly sup- 
presses the chiral condensate. With the proper inclusion 
of the collective baryonic excitations, there is no trace left 
of the chiral first-order transition and the critical end- 
point of the purely mesonic model. The baryon density 
is an order parameter for Nf — N c — 2, and the transi- 
tion line would be expected to give rise to the two-color 
analogue of the baryonic freeze-out. 

The onset of diquark condensation and superfluidity of 
our bosonic baryons, with 50(3) x 50(2) — > 50(3) sym- 
metry breaking at finite quark mass and chemical poten- 
tial, also marks the line at which the the residual 50(3) 
symmetry starts changing in nature from an approximate 
50(5) symmetry as in the normal phase to becoming the 
approximate 50(4) ~ SU(2)l x SU(2)r quasi-restored 
chiral symmetry. Because they are both explicitly broken 
and only approximate symmetries, this vacuum realign- 
ment naturally is a crossover. The quark mass with large 
chiral condensate in the normal phase starts out as a pre- 
dominantly spontaneously generated Dirac mass, and the 
bosonic baryons undergo Bose-Einstein condensation as 
a dilute gas of strongly bound diquarks with the onset 
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FIG. 10. Phase diagram from RG flow with collective bary- 
onic fluctuations (and no chiral I s order transition/critical 
endpoint) compared to the purely mesonic model (A = 0). 

of diquark superfluidity. As their density increases, the 
underlying quark mass rotates into a spontaneous Ma- 
jorana mass leading to a BCS-like pairing. This is the 
relativistic analogue in two-color QCD of the BEC-BCS 
crossover observed in ultracold fermionic quantum gases. 
It is indicated in Fig. [TT] as additional dashed lines in 
the superfluid phase tracing the lines where the quarks' 
Dirac-mass m q = ga equals their chemical potential, i.e. 
fx = m q , see |32j for a comprehensive discussion of this 
crossover within the NJL model. 

In this Figure [TT] we compare the phase diagram of the 
QMD model for two-color QCD as obtained from the full 
RG solution with the mean- field result of Sec. IIIIBI The 
line of the diquark-condensation phase transition, which 
one expects to be of 0(2)-universality, in the QMD model 
RG solution with fluctuations differs more and more from 
that obtained in mean-field QMD and NJL model calcu- 
lations as temperature increases. The first-order transi- 
tion line is washed out by the fluctuations and the as- 
sociated tricritical point as also predicted from next-to- 
leading order x?T [9] turns out to be a mean-field arti- 
fact. As already visible from the T = results for the 
condensates, c.f. Fig. [9j the phase diagrams approach 
one another at small temperatures. 



V. SUMMARY AND OUTLOOK 

In this paper we have developed a Polyakov-loop ex- 
tended quark-meson-diquark model for two-color QCD 
and derived the functional renormalization group equa- 
tion for the grand potential in the leading-order deriva- 
tive expansion. We discussed the mean-field thermody- 
namics of the model and solved the RG flow equation for 
trivial Polyakov-loop, i.e. for the corresponding quark- 
meson-diquark model. In order to correctly describe the 
competing dynamics of collective mesonic and baryonic 



FIG. 11. Comparison of QMD phase diagrams from MF and 
RG calculations, including lines with ga — /i in the superfluid 
phase to indicate the BEC-BCS crossover. 

diquark fluctuations, it is thereby necessary to introduce 
two invariants of the fields in order to account for the rich 
symmetry and symmetry-breaking structure of two-color 
QCD as reviewed in our introduction. The functional RG 
for the effective potential then describes the interplay 
between the collective mesonic and baryonic (diquark) 
fluctuations as summarized once more with showing the 
resulting chiral and diquark condensates over tempera- 
ture and quark chemical potential in a three-dimensional 
plot in Fig. |12| Our numerical solution method on a 
higher-dimensional grid in field space represents impor- 
tant technical progress with many further applications. 



(qq)/(qq>0 
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FIG. 12. Summary: Chiral condensate and diquark conden- 
sate as function of temperature and chemical potential from 
RG calculation. 

One particular advantage of using two instead of the 
usual three colors is that our non-perturbative functional 
methods and model results can be tested against exact 
results and lattice simulations in two-color QCD. First 
important results from such tests are: The expected 0(6) 
scaling at zero density; the relevance of pole masses in the 
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RG framework to correctly describe the onset of diquark 
condensation at the zero temperature quantum phase 
transition of two-color QCD, and the failure of the usual 
screening masses to be capable of that; and finally but 
most importantly, the non-existence of a chiral first-order 
transition and critical endpoint at finite baryon density. 

The latter is not surprising for two-color QCD alone, 
with the BEC-BCS crossover in the superfluid phase of 
the bosonic baryons. We argue, however, that our com- 
parison between the full results with inclusion of collec- 
tive baryonic excitations and the corresponding purely 
mesonic model reveals a general effect, relevant to the 
real world: the chiral condensate drops discontinuously 
at the low temperature liquid-gas transition to nuclear 
matter, and it will continue to decrease with increasing 
baryon density so that one might question whether there 
will be enough chiral-symmetry breaking left for a an- 
other first-order transition at the expected higher densi- 
ties. Similarly, one might speculate that the second-order 
phase-transition line of two-color QCD to diquark super- 
fluidity at finite temperature would lead to the analogue 
of the observed baryonic freeze-out line in the region of 
rapidly increasing baryon density in real QCD. 

Another advantage is that the proper inclusion of bary- 
onic degrees of freedom is much more straightforward 
and much simpler here as in real QCD. While this is 
to a large extent due to the fact that those baryons are 
represented by bosonic diquarks, our study of two-color 
QCD can serve as an important first step towards includ- 
ing diquark-correlations and explicit baryonic degrees of 
freedom in a covariant quark-diquark description by a 
corresponding quark-meson-baryon model for QCD. 

More tests of refined truncations will be performed in 
the future, including a dynamical coupling of the quark- 
meson-diquark model studied here to the full gauge-field 
dynamics of two-color QCD along the lines of what has 
been done already without explicit baryonic contribu- 
tions for QCD [751 UM- 
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Appendix A: Parameter fixing and numerical 
procedure 

In this appendix we briefly outline our parameter fit- 
ting and the numerical methods for the solution of the 
flow equation. 

Since two-color QCD is an unrealistic theory there is 
no canonical choice to fix the parameters to measurable 
quantities. A common approach, often found in the liter- 
ature, is to use the experimentally known N c = 3 values 
and their N c scaling to obtain a consistent N c = 2 pa- 
rameter set. Therefore the assumption f v ~ y/N c yields 
fir = 76 MeV if the usual three-color value f„ =93 
MeV is chosen. Furthermore, we assume that the vac- 
uum pion and sigma masses do not scale with N c and fix 
the pion mass to = m^o = 138 MeV. 

As pointed out above the mass definition at mean-field 
level which is consistent with the Silver Blaze property 
is the pole mass defined via Eq. (34). Usually, in QM 
studies we fix a(T = fi = 0) = f„ together with the 
pion and the sigma masses in the vacuum. The pion and 
sigma pole mass equations Eq. (34) fix the constants A 
and v 2 in the potential V — j — v 2 ) 2 — cu. The ex- 
plicit symmetry breaking constant c is then determined 
by the gap equation. In this way the parameters A, v 2 
and c are found for a fixed momentum cutoff A in the 
vacuum term. As argued above the cutoff A should be 
chosen larger than the largest value of the chemical po- 
tential we are interested in. Here we choose A = 600 
MeV and for comparison A = MeV and adjust 
such that the crossover temperature at fi = coincides 
with the RG calculation. In Table [H] we summarize our 
used parameter values. 



A [MeV] 


m CT [MeV] 


g 


A 


v 2 [10 3 MeV 2 ] 


c [10 s MeV 3 ] 


600 


680 


4.8 


25.05 


-42.710 


2.885 





1055 


4.8 


94.70 


5.575 


1.447 



TABLE II. Parameter values in mean-field approximation 

In the RG setting we adjust the parameter A in the UV 
potential and the explicit symmetry breaking parameter 
c while keeping v 2 — in the UV potential and a Yukawa 
coupling g = 4.8 to match f n and in the IR. As 
there are remaining uncertainties in the determination of 
the pion pole mass via the flow of 2-point function as 
explained in Sec. IV C| we rather determine via the 
onset of diquark condensation at T = 0. 

Finally, we point out some details on the numerical 
procedure to solve the flow equations. The structure of 
the d = flow Q and the full flow for fi = 0, Eq. 



is identical to the flow of the usual three-color QM model 
[B7JI7D]- Several solution methods such as the finite differ- 
ence approach, the Taylor expansion of the effective po- 
tential around a scale dependent minimum or grid tech- 
niques where higher order derivative terms are eliminated 
algebraically are known which all produce consistent re- 
sults. For the full flow equation (49 ) the situation is more 



19 



involved since the effective potential is parametrized by 
two invariants Uk = [/^(p 2 , <i 2 ). 

Here we apply a modified grid algorithm where the 
higher derivatives on the RHS of the flow equation are 
obtained by a two-dimensional-spline fit of the effective 
potential at the respective grid point. The numerical re- 
sults obtained from this procedure agree very well with 
those from an algorithm where the derivatives were ap- 
proximated by finite differences at a fixed discretization 
order. 



Appendix B: Proper-time flow equations 



In this appendix we sketch an alternative derivation 



of the QMD flow equation Eq. ( 49 1 with a proper-time 
regularization, see also 60, 69, 80 and references therein 
for a short introduction and comparison of the proper- 
time with the Wetterich flow. 

Originally, the proper-time renormalization group 
equation (PTRG) was found by an RG improvement of a 
proper-time regularized one-loop effective action |69[ 183] . 
Later, it turned out [5T] that the PTRG flow can be re- 
lated to the Wetterich flow with the background field for- 
malism. The standard PTRG flow can be derived from 
the Wetterich flow when terms proportional to d t T k are 
neglected. The PTRG flow for the QMD model splits 
into a bosonic and fermionic flow dtTk — dfT^B +dtTk.F 
where t = ln(fc/A) denotes the logarithmic RG scale. 
The one-loop expression can be rewritten via Schwinger's 
proper-time representation as 



dr 



dtfairk 2 ) 



(Bl) 



exp (-rrgi) - exp {-tT^)T^ (- M ) 



where the trace runs over momenta and internal indices^ 
As before, the second functional derivative of effective 
action with respect to bosonic (fermionic) fields is de- 

(2) (2) 

noted by T K k ' B (F k F ). The proper-time regulator func- 
tion / a (rfc 2 ) has to fulfill some constrains and the op- 
timal choice, based on incomplete Gamma functions, is 
f a (rk 2 ) = T(a+ l,rfc 2 )/r(a + 1) with a = 3/2 which 
can also be mapped to the optimized regulator in the 
Wetterich flow, for details see [52"] . 



(2) 

In the bosonic case 1^ B is given by Eq. ( 45 1 and can 
be diagonalizcd. Three of the six eigenvalues which are 
related to the three massless pions are degenerate and 



2 For the fermionic part we make use of 75 hermiticity and par- 
ity invarian ce o f the Dirac operator by writing (using the nota- 
tion of Eq. {li} ) A 5 r ( p ) k (cr,n,A; f j,)A 5 = {a, — if, A; -ft) = 

r«W,A;- M ) with^=(^_° 5 ). 



read explicitly p + A 



0J„ 



1,2,3 



where u> n — 2irnT are the bosonic Matsubara frequencies. 
The remaining three eigenvalues p 2 + i = 4, 5, 6 are 
more complicated and are related to the radial er-meson 
and the two diquarks. The three-momentum integration 
separates and can be done analytically. After the proper- 
time integration the bosonic flow is composed of a sum 
over all eigenvalues 



d t T k , B = 

Rewriting 
^ 1 



T k 5 



i=l n£Z 



1 



3 2tt- ^ k 2 + A 



n . k 



i=4 



A « ( W 2 )3 + /3 2 ( w 2 ) 2 +/3 lW 2 +/3 



(B2) 



(B3) 



with the k- and //-dependent coefficient functions a/ and 
Pi listed explicitly in Appendix[D]we arrive at the bosonic 
flow equation Eq. ( [47] ) again. 

In the fermionic sector including the coupling to the 
gauge field via the Polyakov loop variable $ = cos(/3ao) 
we find two 4Nf(— 8)-fold degenerate eigenvalues \ n k of 

the matrix vf F {n)T { ^{-ii) 



X nk = K + ao) 2 -^ 2 +S 2 2 ±2i/VK + a o) 2 + <? 2 |A| 2 , 

(B4) 

with the fermionic Matsubara frequencies v n = (2n + 
1)ttT. This thus yields the fermionic flow of the effective 
potential 



ATk 5 x ^ x - 1 
' ^ ■ A' 

j = ±n£Z n, 



8fc 5 T y. k 2 + g 2 cf> 2 - [i 2 + (u n + a ) 2 

3tt 2 ^ 



n£Z 



(K + ao) 2 + S+ 2 ) ({v n + a ) 2 + E k 2 ) ' 

(B5) 



which reproduces Eq. (48). Evaluating the Matsubara 



sums and combining both contributions then leads to the 
flow equation ( 49 1 . 



Appendix C: RPA Meson/ Diquark polarization 
functions 



For convenience we indicate the explicit expressions 
for the meson/diquark polarization functions for vanish- 
ing spatial external momentum. These can be calcu- 
lated most conveniently using massive energy projectors 
[39] , The Polyakov-loop enhanced quark/antiquark oc- 
cupation numbers N q are defined in Eq. ( 33 1 and reduce 
to the Fermi Dirac distribution for $ = 1. To comply 



with conventions in the literature the polarization func- 
tions are given in a complex basis <p = {?■> ^1 A, A*) and 
correspondingly with Eq. ( 29 ) replaced by 



ny(p) = Tr, 



dT 



(2) 



<^MF 



GMF{p + q)-^j- 



0MF 



Gmf(q) 



(CI) 
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There is a subtlety concerning the u> —> limit at fi- \p\ — >• at lu = 0. This leads to a discontinuity at u> = 

and gives rise to additional contributions (L^n^ for the 
zero mode cj = (see below). Note that polarization 
functions for Goldstone and would-be- Goldstone modes 
are protected from these contributions <5 Wi oII , i.e., only 
n^g. is nonzero in the normal phase. Using m q — ga and 



nite temperature T > as mentioned in Sec. Ill C The 
standard procedure within the imaginary-time formal- 
ism assumes that the external Euclidean pq — — ilj is 
a discrete Matsubara frequency 2imT (n 6 Z). One 
then assumes additional analyticity properties to define a 
unique analytic continuation. The polarization functions 
are then singular in the origin of momentum space and 
one needs to maintain a finite spatial external momen- 
tum p to define dynamic scccning masses via the limit 



<T 



m 2 q , e+ 



6 q ± Ml E q = \/ £q + S^ 2 aS m 

Eqs. ( 18 1 the polarization functions are given bjj^j 



n^(u,T) = -4N c g 2 5, 



E+E- - e+e- - g 2 d 2 /J_ J_ 
u 2 - (E q - E+) 2 \E+ Eg 



13 J (2tt)3 
E+E- + e+e- + g 2 d 2 ( \ 1 \ . 



Nq{Eg) - N q (E+)) 



(C2) 



d 3 <? g 2 



■E+KT - e« e« - 9 2 rf 2 / 1 1 i 

- 3 — 3 Z Hr-~ (^o(^7) " N ZK)) 

0J 2 -(Eq-E+) 2 \E+ Er ,K q q 9 q ' 



E+E' +e+e- +g 2 d 2 ( 1 1 



1 q ' q g 

L0 2 -{Eg +E+) 2 



(jf + ^) (l-N q (E-)-N q (E+)) 



J q ^q 

,2 j2 



2ff^ 1 

;2 - 4Ef 2 El 



(l-2N q (E q t )) 



U AA (u,T) =n A . A .(-w,T) = 47V c5 2 y ^E 



2 Ei± 

Bf+efT^ 1 



^, n° CT (r), 



CJ 2 - AE„ Eq 



(C3) 

5o,, n AA (T), (C4) 



d 

Haa -U:T) : II . A (a.-.D : -4V c5 2 / J 



E 



g 2 rf 2 1 

J 2 - AE+ 2 E± 



(l-2N g (Ef)) 



<L, n AA »(T), 



n CTA (^,T) =n Ao .(w,T) = n aA «(-w,T) = n A . <r (-«,T) 



2e* ± a; 1 
W 2 _ AE f Ef 



(1 - 2W,(£±)) 



■<L,on° A (r), 



with additional contributions for w = of the form, 



]■"• »«2 / pi 2 , ±2 „2j2 n 

d q m q E q +e q - g d 



Et 



n AA (T)=n A . A ,(T) = -2iv c . 



d 3 g 
(2tt) 3 



E 



'K ±2 -e ±2 



E 



±2 



(-2jV^±)) , 
) (-2A^)) , 



3 i 

n Aff (r) =n A , CT (T) = n° A (r) = n° A ,(T) = -2V2N c9 2 [ p>Y,—^h H"W)) , 
n AA ,(r)=n A . A (r) = -27V c5 2 / ^£^(-2^)). 



(C5) 
(C6) 

(C7) 
(C8) 
(C9) 
(CIO) 



3 Note that apart from the missing zero-mode contributions 
5u,onJ., the prefactor of the terms proportional to in n CTCT in 



the last line of Eq. | |C3| l differs by a factor of 2 from the corre- 
sponding terms given in 39 , 40 . 



The additional contributions n (T) vanish for T — >• 
but are required to ensure consistency with the screening 
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mass definition from the effective potential at finite T, 



c.f. Eqs. (361 and (37). 



Appendix D: Coefficients in the bosonic flow 
equation 

In this appendix we list the expressions for the coeffi- 



cient functions on and /3; appearing in Eq. ( 49 1 . 



a 



- 2d 2 U k . dd + Uk. p + 2p 2 Uk, PP ) 



a =3fc 4 + 4fc 2 (-4/i 2 + 2U k ,d + Aa~u k 

+ 4 (4 M 4 + Ul d + 2U k , d {d 2 U kM + U k . p + 2p 2 U k>pp ) - 4/i 2 ([/ M + d 2 U k , dd + U k . p + 
+2d 2 {U kidd U k , p - 2p 2 Ul pd + 2p 2 XJ k . dd U k .pp)) 
1 '- 8U k . d + 8d 2 U k , dd + W Kp + 8p 2 U k , pp 



i =6fc 
a 2 =3 



k,d 



' u k , p 

,,2/ 



2p 2 U k „ 



p 2 =3fc 2 + 8^ 2 + 4U k . d + 4d 2 U k , dd + 2U k>p + Ap 2 U kjPp 



(Dl) 
(D2) 
(D3) 



A) =(k 2 - in 2 + 2U k4 ) (k 4 + 2k 2 {-2p 2 + U kid + -za-u k . dd + t 

+A(-2p 2 U Kp + U kid U ktP + 2d 2 U k4d U k , p - M 2 p 2 Ul pd + 2p 2 (-2p 2 + XJ k4 + 2d 2 U kM )U Kpp )) (D4) 
Pi -3fc 4 + Ak 2 {2U kA + 2d 2 U k4d + U Kp + 2p 2 U k:P p) 

+ 4 (4 M 4 + C/ fe 2 d - Ap 2 (U k4 + d 2 U Kdd - U Kp - 2p 2 U k . P p) + 2U Kd {d 2 U k , dd 
Y2d 2 {XJ kM XJ k ,p - 2p 2 Ul pd + 2p 2 U k . dd U k .pp)) 

, n 2 , 4TT .A llTT , C\T T , A ItT 



(D5) 
(D6) 
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